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ABSTRACT 

The star formation history of nearby early-type galaxies is investigated via 
numerical modelling. Idealized hydrodynamical N-body simulations with a 
star formation prescription are used to study the minor merger process 
(1/10 < M1/M2 < 1/4; Ml < M2) between a giant galaxy (host) and a less 
massive spiral galaxy (satellite) with reasonable assumptions for the ages and 
metallicities of the merger progenitors. We find that the evolution of the star 
formation rate is extended over several dynamical times and shows peaks which 
correspond to pericentre passages of the satellite. The newly formed stars are 
mainly located in the central part of the satellite remnant while the older stars 
of the initial disk are deposited at larger radii in shell-like structures. After 
the final plunge of the satellite, star formation in the central part of the rem- 
nant can continue for several Gyrs depending on the star formation efficiency. 
Although the mass fraction in new stars is small, we find that the half-mass 
radius differs from the half-light radius in the V and H bands. Moreover syn- 
thetic 2D images in J, H, NUV, H/3 and V bands, using the characteristic filters 
of the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST), 
reveal that residual star formation induced by gas-rich minor mergers can be 
clearly observed during and after the final plunge, especially in the NUV band, 
for interacting systems at {z < 0.023) over moderate numbers of orbits 2 
orbits correspond to typical exposure times of ~ 3600 sec). This suggests that 
WFC3 has the potential to resolve these substructures, characterize plausible 
past merger episodes, and give clues to the formation of early- type galaxies. 

Key words: galaxies: formation - galaxies: interactions - galaxies: structure - 
galaxies: kinematics and dynamics - galaxies: photometry - methods: N-body 
simulations 



1 INTRODUCTION 

The formation and evolution of early-type galaxies 
(hereinafter ETGs) remains an outstanding puzzle in 
contemporary astrophysics (see Cimatti 2009 for a re- 
cent review). According to earlier studies, ETGs seem 
to form a class of objects with simple and well defined 
properties: they tend to have a smooth morphology, an 
old stellar population, a red optical colour and are free 
of cold gas and ongoing star formation (Searle, Sargent 
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& Bagnulo 1973; Larson 1975). However, this classic pic- 
ture is challenged by recent observations. For instance, 
the presence of cold gas and dust as well as evidence 
of recent star formation has been detected in several 
systems (see, for instance, Goudfrooij et al. 1994; Mac- 
chetto et al. 1996; Trager et al. 2000; Yi et al. 2005; Mor- 
ganti et al. 2006; Kaviraj et al. 2007; Sarzi et al. 2008, 
Clemens et al. 2009b), and the associated star formation 
rate can even be comparable to that in normal spiral 
galaxies although these events are rather rare (Fukugita 
et al. 2007). Moreover, recent studies based on GALEX 
(Galaxy Evolution Explorer) near-UV photometry of a 
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large sample of ETGs found that 30% of them experi- 
ence ongoing or recent star formation (Schawinski et al. 
2007; Kaviraj et al. 2007), and a similar trend has been 
obtained from mid-infrared Spitzer data in the Coma 
and Virgo Clusters (Bressan et al. 2006; Clemens et 
al. 2009a) or from the COSMOS survey (Kaviraj et al. 
2010). 

Understanding the formation of early-type galaxies, 
and in particular their star formation history, is of cru- 
cial importance for setting strong constraints on models 
of galaxy formation. As has been widely discussed in 
the literature, there are currently two competing sce- 
narios for the formation of ETGs. On the one hand 
there is the monolithic collapse model in which galaxies 
are formed in short, highly efficient starbursts at high 
redshift {z ^ 1) and evolve purely passively thereafter 
(Eggen, Lynden-Beh & Sandage 1962; Larson 1974; Ari- 
moto & Yoshii 1987). On the other hand, the hierarchi- 
cal structure formation model suggests that the most 
massive galaxies are formed from successive mergers of 
smaller entities with accompanying star formation over 
a cosmological time-scale (Toomre 1977; Blumenthal et 
al. 1984; Kauffmann, White, & Guiderdoni 1993). Al- 
though these models are in good agreement with some 
aspects of observational data, they both face some prob- 
lems. For instance, the monolithic model is supported by 
the existence of, and small scatter in, galaxy scaling re- 
lations such as the colour-magnitude relation (Sandage 
& Visvanathan 1978), the fundamental plane (Djorgov- 
ski & Davis 1987; Jorgensen, Franx & Kjaergaard 1996; 
Saglia et al. 1997) and the Mg-cr relation (Colless et al. 
1999; Kuntschner et al. 2001), but cannot easily explain 
the residual star formation and cold gas found in ETGs. 
As far as hierarchical models are concerned, although 
they predict galaxy interaction, they are in disagree- 
ment with the so-called "down sizing" phenomenon. In- 
deed, observations of deep surveys {z > 1 — 2), such 
as the Las Campanas Infrared Survey, HST Deep Field 
North and Gemini Deep Deep Survey (GDDS) have re- 
vealed an excess of massive early-type galaxies under- 
going "top-down" assembly with high inferred specific 
star formation rates relative to predictions of the hier- 
archical scenario (Kodama et al. 2004; Glazebrook et al. 
2004; Cimatti, Daddi & Renzini 2006). 

It is now well known that ETGs have consider- 
able substructure (e.g. from SAURON and GALEX) 
which is interpreted as a result of mergers in the past 
several gigayears. To study this plausible process, we 
have compared the ultra-violet (UV) colours of nearby 
(0.05 < z < 0.06) early-type galaxies with synthetic 
photometry derived from numerical simulations of mi- 
nor mergers, with reasonable assumption for the ages, 
metallicities and dust properties of the merger progeni- 
tors (Kaviraj et al. 2009, hereafter K09). We found that 
the large scatter in the ultra-violet colours of interme- 
diate mass early-type galaxies in the local universe and 
the inferred low-level recent star formation in this ob- 
jects can be reproduced by minor mergers in the stan- 
dard ACDM cosmology. In the present work, our aim is 
to study the evolution of the internal structure of these 
objects using the same methodology but with higher 
mass resolution, in order to help understand in more 



detail the observational signatures of satellite minor 
merger events with different mass ratios, gas-fractions 
and orbital configurations. This work is also motivated 
by the recent installation on the Hubble Space Tele- 
scope (HST) of NASA's Wide Field Camera 3 (WFC3t) 
whose optical design provides a large field of view and 
high sensitivity over a broad wavelength range, excellent 
spatial resolution and a stable and accurate photometric 
performance. It features two independent imaging cam- 
eras, a UV/optical channel (UVIS) and a near-infrared 
channel (IR). Both channels are equipped with a broad 
selection of spectral filters and some of these have been 
considered in our study in order to identify disrupted 
satellites and trace the history of the merger process. 

This paper is organized as follows: in section 2 we 
summarize the numerical methodology that we have de- 
veloped to produce WFC3 synthetic images; in section 
3 we present results on the star formation history of 
nearby ETGs, while the last section presents our con- 
clusions. 



2 NUMERICAL METHODOLOGY 
2.1 Simulations 

The numerical methodology used in the present paper is 
described in detail in K09 where we refer the reader for 
more information. For the sake of clarity, we summarize 
the main steps below. 

Our aim is to study minor mergers between an el- 
liptical galaxy (E) and a satellite (spiral galaxy) us- 
ing idealised N-body+SPH simulations. Both galaxies 
are created following Springel, Di Matteo & Hernquist 
(2005). First, the elliptical is modelled using a spher- 
ical dark matter (DM) halo and a stellar component, 
with a total mass of lO^^M©. Stars contribute 4% of 
this value in agreement with recent observational stud- 
ies (see for instance, Jiang & Kochanek 2007; McGaugh 
et al. 2009). A Hernquist profile reproduces the de Vau- 
couleurs i?^/^ surface brightness profiles of typical el- 
liptical galaxies. The effective radius of the projected 
brightness is re — 4.29 kpc. As far as the satellite is 
concerned, we used a compound galaxy model which 
consists of a spherical DM halo and a rotationally sup- 
ported disk of gas and star (but no bulge) with inde- 
pendent parameters describing each of the structural 
component. The mass of the disk represents 5% of the 
total mass of the satellite and the gas mass fraction in 
the disk is fixed to a constant value 20% (of the total 
mass of the disk), consistent with observed values from 
the SDSS (Kannappan 2004). In all simulations, satel- 
lites are put on prograde or retrograde parabolic orbits 
(Khochfar & Burkert, 2006), with varying pericentric 
distance and initial separations of rini — 100 kpc. 

The simulations are performed using the publicly 
available code GADGET2 (Springel 2005) with added 
prescriptions for cooling, star formation and feedback 
from Type la and II supernovae (SN). It is worth men- 
tioning that gas particles with T < 2 x 10^ number 

t http://www.stsci.edu/hst/wfc3 
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Table 1. Summary of merger simulations. 



run"^ 


iliCtoO LCXljlKJ 


c* 




orbit^ 


Ai 


1:10 


0.01 


8 


prograde 


A2 


1:10 


0.05 


8 


prograde 


A3 


1:10 


0.1 


8 


prograde 


Bi 


1:6 


0.01 


8 


prograde 


B2 


1:6 


0.05 


8 


prograde 


B3 


1:6 


0.1 


8 


prograde 


Ba 


1:6 


0.05 


6 


TirnPTPirlp 


^5 


1:6 


0.05 


4 


prograde 


Bq 


1:6 


0.05 


2 


prograde 


Br 


1:6 


0.05 





prograde 


Bs 


1:6 


0.05 


8 


retrograde 


Ci 


1:4 


0.01 


8 


prograde 


C2 


1:4 


0.05 


8 


prograde 


C3 


1:4 


0.1 


8 


prograde 


Di 


1:1 


0.05 


8 


prograde 



•^Name of run 

^Merger mass ratio (host:satellite) 
"^Star formation efficiency 
"^Pericentric distance of orbit in kpc 
^type of orbit 



A : Projet objects 



Apply to stars 



B : Create spectrum 



Apply to gas 



C : Attenuate spectrum 
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Filter 2 




Filter 3 










E : Map filters to colour 







D : Apply filter 
curves 



F : Export image 



Figure 1. Flow diagram describing the pipeline for the ren- 
dering of telescope views of the objects in the simulation 



density n > 0.1 cm~^, overdensity Apgas > 100 and 
V.v < form stars according to the standard star for- 
mation prescription: dp^/dt = c^pgas/tdyn, where 
refers to the stellar density, tdyn is the local dynami- 
cal timescale of the gas and c* is the star formation 
efficiency. 

In the present study, we use 10 times more particles 
than in K09 in order to more accurately study the evo- 
lution of the spatial distribution of stars. Consequently, 
^ 4, 000, 000 particles are used for each experiment in 
which the particle mass are Mdm = 3.03 X 10^ Mo, 
Mgas = Mdisk = 4.5 X 10^ Mo and Me = 1.35 x 10^ Mo 
for the dark matter particles (DM), gas particles, star 
particles in the disk and star particles in the ellipti- 
cal galaxy (E) respectively. The gravitational soften- 
ing lengths (e) used in each simulation are (physical) 
eoM = = 0.1 kpc and Cgas — eatsk — 0.2 kpc. We 
have also checked that increasing the resolution of the 
simulations as respect to K09 does not affect the star 
formation history for any given experiment which sug- 
gests that the conclusions presented in this work are 
robust. For instance, an example (experiment B2) is ex- 
ibited in Fig. 3 of K09. 

We have considered several realizations whose rel- 
evant free parameters (e.g. initial orbital configuration 
and star formation efficiency) are summarized in Table 
1. Our investigation mainly focuses on minor satellite 
mergers ( 1/10 < mass ratio < 1/4) but we have also 
considered one major merger case (run Di) with which 
to compare the results. 



2.2 WFC3 synthetic images modelling 

This section describes the numerical modelling that en- 
ables the rendering of physically realistic telescope views 
with user- definable filter curves for a simulated instru- 
ment. The rendering process is done using ray tracing 



of light emitted by the simulated objects. For such pur- 
pose, we have developed and used the Alice code"!". 

The instrument view render pipeline is described 
in Fig. 1. Firstly, the star and gas particles are pro- 
jected onto a 2D array corresponding to the instru- 
ment's charge-couple device (CCD) pixels (A). This is 
done using a cosmological distance ladder, assuming a 
ACDM cosmology (e.g. Hogg, 1999), with if = 70 km 
s~^ Mpc~^, Qm = 0.3 and Qa = 0.7 and a flat uni- 
verse (/c = 0). Next, we assign a stellar energy distribu- 
tion (SED) and emission line strength for various emis- 
sion lines to each star particle (B), following Bruzual 
& Chariot (2003). To do this, we use an empirically- 
determined age-met allicity relationship, derived from 
Prochaska et al. (2003). The spectra are frequency- 
shifted according to the Hubble and peculiar velocity 
relative to the simulated observer. Then the code runs 
the emission from the star particles and then the ab- 
sorption from the gas particles. Each of a given pixel's 
spectra and emission line strengths are summed, and 
attenuated (C) according to the gas column density 
(Calzetti et al. 2000, 2001). It is worth mentioning that 
the resolution of the simulations is not enough to resolve 
dense molecular gas or dust. Therefore this may under- 
estimate the total absorption which could affect the pre- 
dicted UV emission. However, we don't think the dust 
attenuation is very high since we showed in K09 that 
the (integrated) UV colours from this set of simulations 
reproduces the observed UV colours of the early-type 
galaxy populations extremely well. If we were system- 
atically underestimating the dust content then we would 
probably see quite a big mismatch between the UV red 
envelope predicted by the model simulations and the 
observed one. In the next step, these spectra are passed 

t www-astro.physics.ox.ac.uk/~samgeen/alice/index.html 
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Table 2. UVIS filters characteristics 



F225W F487N F555W 

0.2 - 0.275 0.48 - 0.49 0.45 - 0.7 /xm 



Table 3. IR filters characteristics 



FllOW F160W 
0.9 - 1.4 1.4 - 1.7 /im 



through a set of transmission filters (D). To do this, we 
linearly interpolate the spectra calculated with the given 
transmission curves. Finally, we produce flux maps cor- 
responding to each transmission fllter and emission line 
(E). We scale the fluxes to the minimum flux detectable 
by the instrument in each wavelength. 

For the purposes of this work, the instrument pa- 
rameters were taken from the WFC3 instrument. For 
instance, we have considered 5 filters: 2 from the near- 
infrared (NIR) field, J (FllOW) and H(F160W), and 3 
from the UV/optical field, NUV(F225W), the narrow- 
band (F487N) and V (F555W) (see tables 2 and 3 
for the associated wavelength bands). The UV colours 
are particularly suitable to detect young stellar popula- 
tions inside interacting galaxies (see for instance, Hib- 
bard et al. 2005) and remains largely unaffected by the 
age-metallicity degeneracy (Worthey 1994) while NIR 
spectral window traces primarily the old stellar popula- 
tions. 

Synthetic 2D images have been produced assuming 
luminosity distances oi D = 100 Mpc and D — 20 Mpc 
corresponding to 2; 0.023 and z ^ 0.005 respectively 
and according the cosmological parameters considered. 
An exposure time, t = 3600 s was used for each of the 
broadband filters, while t = 7200 s was used for the 
narrow-band H/3 filter. We have assumed that each im- 
age is made up of 10 exposures in order to calculate 
an appropriate readout noise. The Field of View (FOV) 
of our synthetic images match those of the WFC3 de- 
tectors, namely: 162 x 162 arcsec for UVIS images and 
133 X 133 arcsec for IR images. All the images have 
been convolved with an appropriate HST Point-Spread 
Function (PSF), derived using the TinyTim^ software 
package (Krist & Hook 2004). Finally, an appropriate 
sky background, Poisson noise and readout noise have 
been added to each of the filtered images. Although 
generated at the full resolution of the UVIS instrument 
(4096 X 4096 pixels) , the UVIS images presented in this 
paper have been re-binned to a resolution of 1024 x 1024. 
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3 RESULTS 

3.1 Evolution of the star formation rates 

In Figs. 2 and 3, we present the time evolution of the 
star formation rates (SFRs) derived from experiments 
described in the previous section. In each case, we have 
let the system evolve until t = 4.0 Gyr from the begin- 
ning of the simulation, which corresponds to ^ 30 tdyn 
where tdyn = {rlir/GMhost(rvir)y^'^ is the dynamical 
timescale at the virial radius of the host elliptical galaxy. 
In our models, rvir ~ 13.5 kpc and Mhost{rvir) ~ 
3.12 X 1O^°M0 which lead to tdyn 1.3 x 10^ yr. This 
value is consistent with observationally determined dy- 
namical times from SDSS parameters (petroRad, pho- 
tometric stellar masses, etc..) that are around 0.1 Gyr 
for ETGs. Young (2002) has also studied interferomet- 
ric images of CO emission in five elliptical galaxies and 
derived some dynamical timescales for gas disks which 
are between 0.07 and 0.2 Gyr. 

We first note the presence of peaks in the evolution 
of the SFRs that are associated with recent pericentric 
passages of the satellite. Indeed, at these specific times, 
an important fraction of available gas is consumed in 
starbursts due to the increase of the gas density during 
the interaction. The last peak in the SFRs is due to the 
final plunge of the satellite remnant. This behavior can 
be clearly seen in Fig. 2, where the evolution of the star 
formation rate derived from experiment B2 is plotted in 
parallel with the variation of the radial distance D of 
the bound part of the satellite remnant (with respect to 
the host galaxy center). We can observe that each SFR 
peak is preceded by a minimum in D. Physically these 
features arise due to tidal forces that induce instabilities 
in the satellite disk, resulting in loss of angular momen- 
tum via dynamical friction and infall to the centre. 

The presence of several peaks in the SFR evolu- 
tion seems to be a general trend in the case of minor 
merger events (mass ratio < 1:4). A similar evolution is 
obtained even while assuming different star formation 
efficiencies or different initial orbital configurations (see 
Fig. 2). We also find similar behaviour for experiments 
B5 and Bq (e.g. for mass ratio 1:6 and rp = 4 and 2 kpc 
respectively). For reasons of clarity, we have omitted 
showing their evolution in Fig. 2. 

However, the result is different for a major merger 
event. The right panel of Fig. 3 indicates that only 2 
peaks are observed. In fact, the timescale of the merg- 
ing process decreases as the mass ratio of the system 
increases. The time delays between the final plunge and 
the first peak are ^ 1.75, ^ 1.25 and ^ 1. Gyr for 
progenitor mass ratios of 1:10, 1:6 and 1:4 respectively. 
This is mainly due to the fact that the satellite galaxy 
gradually loses its energy and angular momentum un- 
der the action of dynamical friction, and finally sinks to 
the center of the host galaxy. This dynamical friction 
operates more slowly for minor merger events (see, for 
instance, Binney & Tremaine 2008). The ages of the in- 
dividual mini-bursts are determined by the orbit and its 
dynamical evolution. 

We also note that if the star formation efficiency 
c* is low, a larger reservoir of gas remains at the time 
of the final plunge. This can sustain star formation in 
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Figure 2. The variations of the star formation rates derived from 1:6 merger mass ratio scenarios. The left panel shows the 
SFRs evolution for a fixed pericentric distance (r^ = 8 kpc) but with different star formation efficiencies c* values. The middle 
panel presents SFRs derived from experiments Bq (retrograde orbit) and (edge-on parabolic merger). In the right panel, we 
have plotted in parallel the SFR derived from experiment B2 and the evolution of the radial distance of the most bound part 
of the satellite. 




t(Gyr) 



Figure 3. The variations of the star formation rates derived from 1:10 (left panel), 1:4 (middle panel) and 1:1 merger mass 
ratio scenarios. 



the remnant for a longer time as is clearly suggested in 
Fig. 2. This result confirms both earlier studies (Larson 
1974) and previous work based on merger simulations 
which show that gas can settle in a self- gravitating cen- 
tral disk which may feed a central black hole (Mihos & 
Hernquist, 1996 ; Barnes, 2002). 



To finish, it should be pointed out that any ini- 
tial transient starbursts and numerous fluctuations seen 
in the star formation rates before the first pericentric 
distance may result from non-equilibrium initial condi- 
tions. Indeed, such starbursts can happen at the begin- 
ning of the simulation since all stars are unaffected by 
feedback and then, most of them may be eligible to form 
stars. The numerous fluctuations seen in SFRs are also 
explained by the fact that a short period is needed for 
type II supernovae to be effective in order to regulate 
the star formation in the galaxy. 



3.2 Evolution of the spatial distribution of 
stars 

In this section, we focus on the evolution of the spatial 
distribution of stars. In our models, stars can be divided 
into 3 distinct groups: those that initially belong to the 
elliptical galaxy (^*), those that were present in the 
disk before the merger (D*), and finally the new stars 
which have formed during the merger process (A^*). We 
are particularly interested in the latter stars. It is worth 
mentioning that we use the same initial age distribution 
as in K09, i.e, we assume that the elliptical is composed 
of old stars of age 9 Gyr whereas the mean age of the 
stellar population in the initial disk is 5 Gyr. The mo- 
tivation of such choices are described in K09. 

Fig. 4 shows the evolution of the projected distri- 
bution (face-on) of all stars in the experiment B2. After 
1.5 Gyr, we notice that the remnant satellite progres- 
sively forms a so-called "shell galaxy", in agreement 
with previous numerical studies (Quinn 1984; Dupraz 
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Figure 4. The evolution of the projected spatial distribution of stars (face-on) from experiment B^. Lines 1 and 2 show all 
stars in grey scale. In the upper panels, the dash blue lines indicate the trajectories of the satellite while the red arrows shows 
the initial spin of the spiral galaxy. Line 3 and 4 represent the stellar ages of the newly formed stars (in Gyr) at the same 
specific times. Each square is 100 x 100 kpc in size. 



& Combes 1986; Hernquist & Quinn 1988; Hernquist & 
Spergel 1992; Hernquist & Mihos 1995). The age of all 
stars is also plotted with a colour coding. At large ra- 
dius, the shell structure is mainly composed of a fossil 
stellar population of accreted stars . Such trends are 
in good agreement with observations of massive ETGs 
(see for instance Clemens et al. 2009b). This is mainly 
due to the dissipationless process and angular momen- 
tum conservation during the merger. The newly formed 
stars are generally formed in the central part of the 
satellite remnant where the density is higher. On the 



contrary, the gaseous component of the satellite falls 
dissipatively into the potential well of the host galaxy 
and, after the final plunge, fuels rapid star formation. 
This latter phase is clearly active in the inner part of the 
system (< 3 kpc), and after 2 Gyr it may be possible to 
measure a gradient in the age of the newly formed stars 
within 40 kpc, as suggested by Fig. 5. Since these 
new stars may have a higher metallicity than the host 
old ones, the presence of a metallicity gradient is also 
expected, in agreement with observations (Clemens et 
al. 2009b). 



© 0000 RAS, MNRAS 000, 000-000 



Composite star formation histories of ETGs 7 



4.0 
3.5 
^ 3.0 
S 2.5 

D) 

2 2.0 
CO 1.5 
1.0 



0.5 



J 1 M 1 


MM 


MM 


MM 


M M M M U 










■"^=4.0 Gyr= 










t 2 5 Gyrz 










^ t=2.0 Gyr- 










t=1.5 Gyr3 


Sum 


MM 


MM 


MM 


M M M M F 



10 20 30 40 

R (kpc) 



50 60 



Figure 5. Evolution of the distribution of ages of newly 
formed stars (face-on) at different times from the beginning 
of the simulation for experiment B2. 



3.3 Size-Evolution of wet merger remnants 

Simple arguments based on the virial theorem allow us 
to calculate the size evolution of the collisionless compo- 
nent of the merger remnant (e.g. Bezanson et al. 2009; 
Naab et al. 2009). Assuming the individual stellar com- 
ponents of the galaxies are in virial equilibrium: 



fGMf 
2 n 



(1) 



The factor / depends on the individual functional form 
of the profile and is close to / = 0.5 in the case of de 
Vaucouleurs and exponential profiles. The radius r is the 
half-mass radius of the galaxy and M its stellar mass. 
The merging galaxies are set up to be on parabolic or- 
bits with Eorb = 0, which allows us to write down the 
internal energy of the remnant assuming energy conser- 
vation as 



El + E2 + Eor 



io 



Mt ^ Mj 
ri r2 



(2) 



Assuming the remnant will settle down in virial equilib- 
rium, then basic manipulations of Eq. 2 lead to 



(Ml + M2)' 



Mt 
ri 



+ 



Ml 
r2 



(3) 



In the following, we investigate the merger remnants 
and characterize the deviations from this simple virial 
argument based on the rate of dissipation and star for- 
mation occurring during the merger. For this purpose, 
we focus again on experiment B2 which can be consid- 
ered as a typical minor merger scenario. In this case, we 
obtain rrem ^ 5.67 kpc. 

From the simulation, the half-mass radius Re at a 
given time is defined to be the radius of the circle con- 
taining half of the projected stellar mass. We have also 
studied the evolution of effective radius RY , and R?^ 



from V and H bands respectively. RY and R^ are de- 
fined to be the radius containing half of the projected 
flux in each specific band. In order to take into account 
the effects of viewing angle, we have considered 100 dif- 
ferent lines of sight chosen randomly and derived the 
mean value of Re, RY and R^ . 

Fig. 6 shows the evolution of Re, RY , and R^ after 
the final plunge for experiment B2 . The evolution of Re 
tends to a constant value {Re 5.05 kpc) which is lower 
than the predicted value rrem 5.67 kpc. This result is 
in agreement with the study of Covington et al. (2007) 
who show that the model systematically over-predicts 
the radii of remnants and using a wide range of merger 
simulations which account for dissipative energy losses 
and star formation. 

The evolution of RY is different and presents two 
phases. Firstly, a rapid decrease just after the final 
plunge (t ^ 1.6 Gyr) which continues for ^ 0.5 Gyr. 
This period corresponds to a phase of high star forma- 
tion induced by the final plunge as shown in Fig. 2. The 
newly formed stars, mainly located in the central re- 
gion, have a significant contribution of the total flux in 
the V-band leading to the reduction in RY . In the sec- 
ond phase, beginning at t ^ 2.25 Gyr, RY increases due 
to the fact that the stellar population becomes progres- 
sively older. Indeed, such behavior can be clearly seen in 
the lower panel of Fig. 6 where we have plotted the evo- 
lution of the half-mass radius for each stellar component 
E:^ , and . The half-mass radius of the host ellipti- 
cal galaxy tends to a constant value ('-^ 4.2 kpc) which 
is actually very close but lower than the initial value 
(4.29 kpc) and probably due to the adiabatic contrac- 
tion induced by the infall of the satellite at the center. 
The half-mass radius of stars of the initial disk also 
tends to a constant value ( ^ 11.5 kpc), which is quite 
high but explained again by the dissipationless process 
and angular momentum conservation during the merger. 
Finally, the half-mass radius for the newly formed stars 

is rapidly decreasing until t ^ 3.0 Gyr due to the 
high star formation induced after the final plunge. After 
this time, the increase observed in the evolution of RY 
results from the ageing of the stellar population . 

The evolution of R^ follows the same trend, though 
the decrease after the final plunge is less pronounced but 
continues for a longer time (~ 2.0 Gyr). The value of 
R^ is close to Re values but still smaller. It is expected 
that both RY and R^ asymptotically tend to the value 
of Re at late time. We infer that because of even modest 
amounts of star formation, the half-mass radius differs 
significantly from the half-light radius especially in the 
V-band. The observed half-light radius is expected to 
be smaller than the half-mass radius. 



3.4 Evolution of WFC3 bands 

In order to produce synthetic images, we have assumed 
our galaxy system to be in the local universe. Such a 
choice is motivated by the fact that up to two-thirds 
of nearby ETGs contains shells and ripples (see for in- 
stance Malin & Carter 1983; van Dokkum 2005) and 
these fine structures and disturbed morphologies often 
coincide with signatures of recent star formation (e.g. 
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Figure 6. The upper panel shows the evolution of the half- 
mass radius (green line) and effective radius from V (blue) 
and H (red) bands for experiment B2. In the lower panel, we 
show the evolution of the half-mass radius for each stellar 
component: stars that initially belong to the host elliptical 
galaxy (£^*, red line), those that were present in the disk 
before the merger (D*, green line) and the new stars which 
have been formed during the merger process (A^*, blue line). 
In each panel, the horizontal dashed line is the predicted half- 
mass radius value after the merger from arguments based on 
the virial theorem whereas the vertical dashed line is the 
epoch of the final plunge. 



Schweizer et al. 1990; Schweizer & Seltzer 1992). It is 
then particularly interesting to try to characterize this 
spatially resolved fine structure whose detection is possi- 
ble only in our local neighbourhood. For instance, using 
N-body simulations, Feldmann, Mayer & Carollo (2008) 
have shown that this broad tidal features cannot be re- 
produced by equal-mass dissipationless merger but is 
well explained by the accretion of disk-dominated galax- 



ies. Moreover, minor mergers are expected to be more 
common than major mergers at low redshifts (Khochfar 
& Silk 2006; Genel et al. 2008; Khochfar & Silk 2009) 
and such statements support our choice of merger mass 
ratio values. 

We start our investigation with experiment B2 . The 
galaxy pair is supposed to be at a redshift z ^ 0.023 (or 
equivalently at a luminosity distance of D = 100 Mpc) 
in order to facilitate comparisons with future observa- 
tional data. In the following, the magnitude units given 
as measures of spectral flux are W.m~^, and in all syn- 
thetic images, we use a logarithmic scale. 

In Fig. 7, we present the grid map derived by our 
numerical modelling through the J and H bands. It ap- 
pears that it is possible to resolve the host elliptical 
galaxy and the satellite remnant, in particular the shell 
structure which is mainly composed of old stars. How- 
ever, this latter tends to disappear at t — 4.0 Gyr. We 
also note that there are no significant differences be- 
tween J and H bands. Fig. 8 shows the evolution of de- 
rived synthetic images from NUV, H/3 and V bands. The 
ongoing star formation regions can be clearly followed 
through the three filters. In particular, after the final 
plunge, ongoing star formation located at the center of 
the galaxy can be clearly observed. In fact, when com- 
paring with Fig. 4, only stars with ages < 0.5 Gyr can 
be identified. Thus, the combination of IR and UVIS im- 
ages allows us to separate different stellar populations 
and then distinguish the most bound part of the satellite 
remnant, composed of young stars, from the host galaxy, 
composed of older stars. This combination also gives 
useful clues on the formation of ETGs: while the shell 
structure revealed in IR images support a past merger 
scenario, evidence or not of the presence of young stars 
in UVIS images bring additional constraints on the wet- 
ness/dry ness of the merger. 

The same study has been performed for objects as- 
sumed to be at 2; 0.005 {D = 20 Mpc). In this case, 
only the inner part is imaged due to the instrument 
FOV, as shown in Figs 9 and 10. The evidence of recent 
star formation is again clearly visible especially in the 
NUV band after the final plunge (t > 1.5 Gyr). Some 
structures can also be distinguished around the central 
part and represent signatures of a recent merger event. 

We have also studied how the free parameters (star 
formation efficiency and initial mass ratio) affect our re- 
sults. For instance. Fig. 11 shows the evolution of syn- 
thetic images from the J and NUV bands in the experi- 
ment Bi. In this case, although the star formation effi- 
ciency is lower, the satellite remnant can also be identi- 
fied in the NUV band. In the J band, the shell structure 
is less pronounced. 

Comparison with mass ratios 1:4 and 1:10 presents 
similar conclusions although the structures are better 
resolved in the case of the 1:4 mass ratio due to the 
higher star formation induced during the merger (see 
Fig. 12 and 13). for instance maximum values in the 
NUV band between 1:4 and 1:10 merger mass ratios 
differ by up to an order of magnitude while they are 
comparable in the H band. However, simulations with 
the 1:4 merger mass ratio should be regarded as a lim- 
iting case, since it is close to the major merger scenario. 
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Figure 7. The evolution of synthetic images from experiment through J (first hue) and H band (second hue) assuming the 
observed system at z ~ 0.023. The images have been Gaussian smoothed with a sigma of 2 pixels. 
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Figure 8. The evolution of synthetic images from experiment B2 through NUV (first line), (second line) and V (third line) 
band assuming the observed system at z ~ 0.023. 



4 DISCUSSION 

The formation of early- type galaxies remains a major 
unresolved problem in astrophysics, and in particular 
their star formation history is known to be more com- 
plex than that of a single population model. Studies of 
nearby galaxy populations can provide useful clues to re- 
solving this problem, by setting tight constraints on the 
level of the most recent star formation activity. In the 
present paper, we have used idealized hydrodynamical 



simulations, including a standard star formation pre- 
scription, to study the minor merger process between 
an elliptical and a spiral galaxy which is expected to 
be the dominant channel in terms of the morphology of 
the merging galaxies (Khochfar & Burkert 2003). There 
is now observational evidence that massive ETGs may 
have undergone not more than one major merger since 
z - 0.7 (Bell et al. 2006; Mcintosh et al. 2008). Such re- 
sults are supported by theoretical evidence that massive 
halos (> 10^^ M0) experience only one major merger be- 
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Figure 9. The evolution of synthetic images from experiment B2 through J (first hue) and H band (second hue) assuming the 
observed system at z ~ 0.005. 
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Figure 10. The evolution of synthetic images from experiment B2 through NUV (first line), (second line) and V (third 
line) band assuming the observed system at 2; ~ 0.005. 



tween z = 2 and 2; = 0, while minor mergers are much 
more common (Khochfar & Silk 2006; Genel et al. 2008; 
Khochfar & Silk 2009). We use idealised simulations in 
this work for two main reasons. Firstly, it allows us to 
reach a high mass resolution lO^^^M©) with rea- 
sonable computational costs compared to cosmological 
simulations. Secondly, it permits us to consider a wide 
range of free parameters (including star formation ef- 
ficiencies and initial orbital parameters) and to study 
how the results change accordingly. 



We have generated 15 experiments by considering 
different merger mass ratios, (1:10, 1:6, 1:4 and 1:1) 
and three star formation efficiencies (c* =0.01, 0.05 and 
0.1). We find that the satellite is progressively disrupted 
by tidal torques and the star formation induced by the 
merger process shows several peaks in its evolution, as- 
sociated with the passage of the remnant at decreas- 
ing pericentric distances. The satellite remnant progres- 
sively forms a shell structure within the host galaxy, 
which is composed of old stars at large radii and is ex- 
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Figure 11. The evolution of synthetic images from experiment Bi (mass ratio 1:6 and c*=0.01) through J (first hue) and NUV 
(second hue) assuming the observed system at 2; ~ 0.023. 



1 1 1 r] 1 1 M 1 M 1 1 

-19.2 -18.8 -18.5 -18.2 -17.9 

J' 

0.5 Gyr J (FllOW) 


1 1 1 ^ 1 1 M 1 1 M 1 

-19.2 -18.8 -18.5 -18.2 -17.9 

1.0 Gyr J (FllOW) 


1 1 1 n 1 1 M 1 M 1 1 
-19.2 -18.8 -18.5 -18.2 -17.9 



1.5 Gyr J (FllOW) 


1 1 1 1 1 M 1 M 1 1 
-19.2 -18.8 -18.5 -18.2 -17.9 

o 

2.5 Gyr J (FllOW) 


1 1 1 [' 1 1 1 M 1 M 1 1 
-19.2 -18.8 -18.5 -18.2 -17.9 

4.0 Gyr J (FllOW) 


1 1 1 [ 1 1 1 1 1 1 1 1 1 1 

-22.9 -21.9 -20.9 -20.0 -19.(1 

0.5 Gyr NUV (F225W) 


1 1 1 1 ] 1 1 1 1 1 1 1 1 1 

-22.9 -21.9 -20.9 -20.0 -19.« 

1.0 Gyr NUV (F225W) 


-22.9 -21.9 -20.9 -20.0 -19.(1 

1.5 Gyr NUV (F225W) 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^i^^ 

-22.9 -21.9 -20.9 -20.0 -19.0 

2.5 Gyr NUV (F225W) 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 pM^^ 

-22.9 -21.9 -20.9 -20.0 -19.0 

4.0 Gyr NUV (F225W) 



Figure 12. The evolution of synthetic images from experiment A2 (mass ratio 1:10 and c*=0.05) through J (first line) and 
NUV (second line) assuming the observed system at z ~ 0.023. 



plained by the fact that the stellar accretion and en- 
suing relaxation occurs via dissipationless mechanisms. 
On the contrary, the center is dominated by the newly 
formed stars due to the infall of the gas to the cen- 
tre in a dissipative process. A positive radial age gradi- 
ent is clearly seen and remains constant in time. All of 
these results are consistent with the conclusions of ob- 
servational studies of massive ETGs (see for instance, 
Clemens et al. 2009b). 

These typical shell structures are of particular in- 
terest as a means of characterizing the last merger expe- 
rienced by the host galaxy, i.e, the mass ratio, dryness 
or orbital configuration. For instance, recently Canalizo 
et al. (2007) have tried to set constraints on the age 
of the plausible merger responsible for shell structures 
observed in a quasar host galaxy from HST/ACS. By 
using N-body simulations, they found that the merger 
occurred a few hundred Myr to ^ 2 Gyr ago when 



matching the observed and simulated outermost shell 
and causal connection with the quasar activity. But the 
accurate estimation of the age of the merger proves to 
be a very hard task due to the high number of orbital 
parameters, even when combining dynamical and photo- 
metric properties of the galaxy. For instance, Di Matteo 
et al. (2009) investigated the evolution of star forma- 
tion efficiency in galaxy interactions and mergers using 
numerical simulations and found that retrograde orbits 
or high pericentric distances favor the enhancement of 
the star formation. Such degeneracy may be difficult to 
deal with. 

Our main result is to demonstrate that colour infor- 
mation provides an important additional constraint via 
the induced star formation that together with the asso- 
ciated morphology provides further and important in- 
sights into the merger parameters. In particular, future 
observations especially with the resolution and NUV 
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Figure 13. The evolution of synthetic images from experiment C2 (mass ratio 1:4 and c*=0.05) through J (first fine) and NUV 
(second fine) assuming the observed system at z ~ 0.023. 



power of WFC3 will open new vistas on star formation 
histories of ETGs. We have derived realistic synthetic 
images through J, H, NUV, and V filters that can be 
compared with future WFC3 observations. Comparison 
of synthetic maps in the NUV and J bands for differ- 
ent merger mass ratios and star formation efficiencies 
gives similar results: evidence of recent star formation 
can be clearly detected before and after the final plunge, 
especially in the NUV band. 

Our study does not include AGN feedback. We will 
consider this in a future paper. Here we note that while 
such a mechanism is supposed to play an important role 
in the early phases of the spheroid star formation, its 
role in the late phase is still poorly known. As suggested 
by our results, the infall of the gas at the final plunge can 
settle in a self- gravitating central disk which may feed a 
central black hole. For instance, N-body simulations in- 
cluding feedback from black hole accretion indicate that 
the quasar activity is coincident with the coalescence of 
the nuclei and eliminates all star formation afterwards 
(Springel, Di Matteo & Hernquist 2005). This will prob- 
ably lead to a reduction of the total flux in the UVIS 
bands. 

In summary, minor mergers seem to play a funda- 
mental role in the late phases of ETG evolution. Minors 
mergers may be the principal mechanism behind the 
large UV scatter and associated low-level recent star for- 
mation observed in ETGs in the nearby universe (Kavi- 
raj et al. 2009). In general, we find composite popula- 
tions contributing to the star formation rate which per- 
sists for 2-3 Gyr, corresponding to several orbital times. 
Minor mergers naturally lead to star formation that is 
extended over a much longer period than for the single 
population starburst. The limiting case of major merg- 
ers results in much shorter star formation time-scales. 
Size evolution is another natural consequence of the mi- 
nor merger model. The considerable amount of dissi- 
pation and central star formation reduces the apparent 
half-light radius of the merger remnant compared to its 
half-mass radius (see Fig. 6), thus indicating that size 



determinations based on rest-band magnitudes sensitive 
to star formation will clearly underestimate sizes by fac- 
tors of 1.2-1.6 for a few Gyr after the merger. Taking our 
results at face value and considering that the frequency 
of gas-rich minor mergers is larger at earlier times, the 
size-evolution of elliptical galaxies from high to low red- 
shift could be milder than assumed so far. 

To conclude, the present work shows that minor 
mergers induce amounts of star formation in ETGs 
which can be measured through UV bands. WFC3 
represents the best instrument to study these minor 
merger events because it has a matching UV and op- 
tical FOV and gives the resolution to see young sub- 
structures (which would not be possible with GALEX 
for instance). The previous ACS/HRC UV detectors had 
a tiny FOV so it was not possible to study a galaxy up 
to 1 effective radius because one would get the whole 
galaxy in the optical image and only a fraction of its 
core using the UV. With WFC3, it is now possible to 
map the entire galaxy in both the UV and optical mak- 
ing it possible for the first time to perform spatially re- 
solved star formation histories in ETGs at low redshift 
using the UV. Moreover the ability to see the young 
substructures is important because it rules out UV flux 
from old sources such as horizontal branch stars (which 
would follow the optical light profile). 
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